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\Q • In the present work, inhomogeneous Tolman-Bondi type dust space-time is studied on the brane. 

There are two sets of solutions of the above model. The first solution represents either a collapsing 
model starting from an infinite volume at infinite past to the singularity or a model starting from a 
singularity and expanding for ever having a transition from decelerating phase to accelerating phase. 
The first solution shows that the end state of collapse may be black hole or a naked singularity 
depending signs of various parameters involved. The second solution represents a bouncing model 

Fj^ ■ where the bounce occurs at different comoving radii at different epochs. 

t^ - ' PACS numbers: 04.20 -q, 04.40 Dg, 97.10. CV. 
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In the seminal papers [1, 2] the introduction of branes into cosmology offered a novel approach to the 
understanding of the evolution of our universe. It was proposed that our 4D universe is said to be a singular 
hypersurface, a 3-brane embedded in a (4 + d) dimensional space-time, which is called the bulk. Here in 
the bulk, we have considered no matter field. The matter field is all confined to the brane and the gravity 
propagates in the extra dimensions as well. The effective equations for gravity in four dimensions were obtained 
by Shiromizu et al [3] and discussed by Roy Maartens [4]. Here Israel's boundary conditions were used along 
, with Z2 symmetry for the five dimensional bulk space-time embedding our brane universe at a fixed point 
y = yo. If we add a cosmological constant in the bulk, the solutions of the Einstein's equations G a b = n\T a i, can 
be obtained, where the universe starts with a non-conventional phase and then enters the standard cosmological 
phase [5 - 11]. Friedmann like equations are usually solved on the brane for matter localized on the brane [12]. 
The usual case in the brane model recollapses if it is embedded in an anti de-Sitter (Ag < 0) five dimensional 
bJ)' bulk, where the extra dimension is space-like. In a recent work of Ponce de Leon [13] the acceleration of the 
present universe is studied as a consequence of the time evolution of the vacuum energy in cosmological models 
• i-h based on brane world theories in 5 dimensions. The basic two assumptions in this paper are that the universe 
is spatially flat and the matter density decrease as an inverse power of the scale factor. If A5 = then the 
induced cosmological term in 4D is positive and the brane universe expands continuously in a power law time 
dependence. Again for A5 > 0, this universe becomes dominated by a positive cosmological term A4. The effect 
is an asymptotic de-Sitter expansion which occurs regardless of the signature of the extra dimension. On the 
otherhand, with a single timelike extra dimension the cosmological brane world model experiences a natural 
bounce without ever-reaching the singular state [14]. 



In the present work, we have considered the Tolman-Bondi space-time on the brane. This is an inhomoge- 
neous and anisotropic space-time and differs from the homogeneous Friedmann space-time so far mentioned 
above regarding their dynamical behaviour. It is perhaps for the first time that Tolman-Bondi type dust 
space-time is studied on the brane. We'll consider only the spatially flat model (f(r) = 0), which is consistent 
with the present observational data. The projection of the bulk Weyl tensor on the brane is also assumed to 
be zero. There occur two different cases so long as we assume the anti de-Sitter bulk and the extra dimension 
to be space-like. In the case I, the brane world collapses from an infinite past to a singularity at a finite epoch 
or otherwise begins from the singularity and expands indefinitely with deceleration near the big bang and 
acceleration at the late stage. In case II, the brane world model has a lower bound and hence singularity-free, 
similar to the case mentioned above (Ref. [14]). Here we get a contracting model bouncing from a lower bound, 
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which however takes place at different shells with different comoving radii at different epochs. This is usual for 
an inhomogeneous space-time. It is to be noted further that in the Tolman-Bondi space-time on brane under 
consideration, there is no recollapsing model. 

Field equations in the Brane: 

The modified 4D Einstein's equations on the brane embedded in the 5D bulk are usually presented in the 
form [1, 3] 



Gfiv — —h-XQuv + KjjT^v + K^S^v — (1) 

where A4 is the effective four dimensional cosmological constant, K5 and K4 are the coupling constants related 
to the gravitational constants in 5D bulk and 4D brane respectively. They arc further related with the funda- 
mental 5D Planck mass and the effective Planck mass on the brane. T Mi , is the usual matter energy tensor on 
the brane. The two additional terms and E^ v on the r.h.s of (1) stand respectively for the local quadratic 
energy-momentum correction term and the free energy induced on the brane by the 5D Weyl tensor in the bulk. 
The local correction term arises from the extrinsic curvature of the brane. Their explicit expressions are given 

by 

V = T M „ - \T«T ua + l^(3T^T Q/3 - T 2 ) (2) 
Finally, E^ u in (1) is the projection of the bulk Weyl tensor on the brane defined by 



E ab = C abcd n c n d (3) 

Here C a bcd is the usual Weyl tensor in the bulk, n a is the unit normal to the hypersurface y = constant, 
representing the brane world, where y represents the extra fifth dimension. In the above the latin indices refer 
to the 5D bulk and greek indices correspond to the brane world. The four dimensional A4 and 5D cosmological 
constant A 5 [2, 3] are related by the following relation 



K 2 K 2 \ 2 

A 4 = f(A 5 + ^) (4) 

where A is the so-called brane tension or the vacuum energy of the brane world. This brane tension is further 
related to the 4D gravitational constant by a relation of the form 



In this work, we assume that our brane model embedded at y = constant in the conformally flat bulk (zero 
Weyl tensor) is in homogeneous and an isotropic. As a result the term E^ is neglected and is justified in view 
of the late time dominance of the vacuum energy compared to the brane matter density p i.e., A >> p [4] being 
compatible with observations. Using spherically symmetric form the metric is given by 



ds 2 = dt 2 - e a dr 2 - R 2 dn 2 (6) 

where a — a(r, t) and R — R(r, t). 

Now considering matter in the form of dust on brane the independent field equations are given by 

1 a ( n R" R' 2 / R'\ R 2 ■ R » p 2 
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1 n, R'^ „ R R^ P 2 /„N 

R 2 ~* V +2 R + R*= A -k ® 



,'R" a' R'\ a a 2 aR R , p 2 

-^ IT-TIT +7T + T+oo+o= A -^ (9) 



and 



R 2R 2 A 2R R 2A 



R' . R' 



2^-«^=0 (10) 



It follows immediately from the equation (10) 

R n 



1 + f(r) 



(11) 



In general there may be three different choices f(r) — 0, j(r) > 0, f(r) < 0. For simplicity f(r) is assumed 
to be zero (marginally bound case). In this particular case the field equations (7) to (9) simplify to 

R 2 _ R R! p 2 

R 2+2 RW= A + P+ k (12) 



„ R R 2 O 2 r~„\ 

2 R + R 2 = A -TX ™ 

R' R R 1 R p , . 

Rf + RRf + R =A ~2X (U) 

It is to be noted that the energy density enters the field equations quadratically in contrast with the usual 
general relativistic 4D equations. 

Subtracting (13) from (14) we get 

R! R R! R R 2 . . 

K + R~R~R~R 2 ~ = °> (15) 



which can be integrated twice to yield 



R 2 = /ii? 2 + W- (16) 
R 

In (16), /i is in general a function of time and D(r) is a function of radial co-ordinate alone. In order to 
integrate equation (16) we assume further (i = constant, which in fact has two groups of solutions. Also from 
the field equations (12) - (14) using (16), we finally get 

D 2 (r) 

P = 6 ^ 2A + W (17) 



and 

-,2 



| X =A-3/x. (18) 
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Solutions: 

There are three possible types of solutions of the evolution equation (16) depending on the signs of the 
parameter \x and the function D(r) namely, (i) \i > 0, D(r) > 0, (ii) fi > 0, D(r) < 0, (Hi) \x < 0, D(r) > 0. 
However, the third possibility will not be discussed as it gives oscillatory solution, not of much interest in the 
present context. 

Case I : /j, > 0, D(r) > : 

The solution for the scale factor R from equation (16) gives 

R = R (r) Sinh 2 ' z u (19) 
where, Rq — ^jp-, u = a(t a (r) —t), a = ^[Jx and t (r) is an arbitrary function of the radial co-ordinate r. 
Differentiating R we get in turn 



2aR Coshu , - 2a 2 R \Cosh2u - 21 
R = — — — — and R = 



(20) 



Here there may be two situations. The first case corresponds to t < to. The collapse starts from the infinite 
past (t — > — oo) i.e., when R — > oo and finally converges to the singularity R — > at the time to(r'). The rate of 
contraction at both ends is infinitely large. It is a collapsing model of the brane. One should note however that 
different shells collapse at different epochs to(r) because of inhomogeneity. The time of formation of central 
singularity is given by [15 - 17] 



t s = ti H — H m gfafo i 
a r -> 

where we have chosen R — r at t — ti as an initial condition. 



fir 

W) 



1/2 



(21) 



Now at the time of formation of apparent horizon (t a h(r)) for a shell having comoving radial co-ordinate r, 
we must have R 2 = 1. Using dynamical equation (16) it results 

fiR 3 - R + D(r) = 0. (22) 
Thus the time difference between the formation of apparent horizon and the central singularity is [15 - 17] 



tah(r) - t s 



1 



Sink- 1 ( —J -Sink- 1 



Rah 

Ro 



3/2' 



1 lim n Sinh-i(-^ 
a r->0 ^£)( r ) 



1/2 



(23) 



where R a h is the positive real root of the above cubic equation (22) in R. Further, in order to have the 
collapsing process starts from a regular initial data one should have the density to be finite at r = when 
t = ti . Thus near r = the function D(r) should have the power series expansion 



D(r) = D r 3 + D x r A + D 2 r 5 
Then the above time difference simplifies to 



(24) 



tah(r) - t s 



2LWM + D 



r + 



(2fj, + 3£>q)J^ - 4D Q D 2 (fi + D ) 
8D 2 (^ + D r/ 2 



0(r 3 ) 



(25) 



Thus we note that the sign of this time difference depends on various parameters involved. Therefore, 
the spherically symmetric space-time shows collapse to a black hole or a naked singularity near the centre 
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Fig. 1 shows the behaviour of R of eq.(20) with time in case I for fj, = 1, Do = 1 and to = 1. 




t 

Fig.2 

Fig. 2 shows the behaviour of R of eq.(26) with time in case II for fj, = 1, i?o = 1 and to = 1- 

depending on the nature and properties of certain parameters. 

On the otherhand, if t > to, R(t) starts from the singularity R — at t — to and expands indefinitely. 
Although this shows a monotonic expansion R < initially and subsequently changes sign (see figure 1) to 
R > at the later stage. Thus there is a transition from deceleration to acceleration during the evolution of 
the brane model. This result seems to be interesting because the model is consistent with the possibility of the 
structure formation in the early stage (R < 0) and the accelerated expansion in the later stage getting support 
from the observational data. 

Case II : \i > 0, D(r) < 0: 

In view of (16) this case is restricted by the condition that fiR? > \D(r)\. The solutions are given by 



() 



R = R {r)Cosh 2 ^u (26) 

Hence, 

■ 2aR Q Sinhu ~ 2a 2 R [Cosh2u + 2} 

R ZCosh^u ^ R 9Cosh^u (27) 

In this case, R is never zero, which shows that the brane world model is singularity free. This conclusion 
follows immediately from the positive value of R throughout the evolution and hence there is a minimum only 
and no maximum at any stage. R has an infinite magnitude at t — > — oo, subsequently contracts to R = Rq 
when t = to and then again explodes to infinity at t — > oo (see figure 2). So this particular brane world model 
shows bounce. However, it should be noted that due to inhomogeneity in the energy distribution the bounce 
occurs at different comoving radii at different epochs. 
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